Theoretical modeling of spatial and temperature dependent exciton energy in coupled 
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Motivated by a recent experiment of spatial and temperature dependent average exciton energy 
distribution in coupled quantum wells [S. Yang et al., Phys. Rev. B 75, 033311 (2007)], we 
investigate the nature of the interactions in indirect excitons. Based on the uncertainty principle, 
along with a temperature and energy dependent distribution which includes both population and 
recombination effects, we show that the interplay between an attractive two-body interaction and 
a repulsive three-body interaction can lead to a natural and good account for the nonmonotonic 
temperature dependence of the average exciton energy. Moreover, exciton energy maxima are shown 
to locate at the brightest regions, in agreement with the recent experiments. Our results provide 
an alternative way for understanding the underlying physics of the exciton dynamics in coupled 
quantum wells. 

PACS numbers: 71.35.Lk, 71.35.-y, 73.20. Mf, 73.21.Fg. 



I. INTRODUCTION 

An exciton is a bound pair of an electron and a hole 
in a semiconductor. At low densities, exciton are Bose 
quasiparticles with a small mass, somewhat similar to 
hydrogen atoms. At low temperatures, of the order 
of 1 Kelvin, the excitons are expected to realize the 
phenomenon of Bose-Einstein condensation (BEC).^ Re- 
cently, it has been shown that indirect exciton (spatially 
separated electron-hole pairs) in coupled quantum wells 
(CQW) can have a long lifetime and high cooling rate. 
With these merits, Butov et al. have successfully cooled 
the trapped excitons to the order of 1 Ki^ Although there 
is not enough evidence to prove that these excitons are 
condensed into the BEG state, it is fascinating enough to 
observe several novel features of the photoluminescence 
(PL) patternsM 

Major features of the macroscopically ordered exciton 
states observed are summarized as follows.— (i) Two ex- 
citon rings are formed. When the focused laser is used to 
excite the sample and prompt luminescence is measured 
in the vicinity of the laser spot, a ring, called the internal 
ring, is formed. A second ring of PL appears away from 
the source at the distance about 1 mm, called the exter- 
nal ring, (ii) The intervening region between the internal 
and external rings are almost dark except for some local- 
ized bright spots, (iii) Periodic bright spots appear in the 
external ring. The bright spots follow the external ring 
either when the excitation spot is moved over the sam- 
ple, or when the ring radius is varied with the excited 
power, (iv) PL is washed out eventually with increasing 
temperatures, (v) In an impurity potential well, the PL 
pattern becomes much more compact than a Gaussian 
with a central intensity dip, exhibiting an annular shape 



with a darker central region, (vi) With the increase of the 
laser power, exciton cloud first contracts then expands, 
(vii) Exciting sample by higher-energy lasers, the dip can 
turn into a tip at the center of the annular cloud. 

The above experimental facts raised several interest- 
ing questions which need to be clarified, (i) What is the 
reason behind the two-ring structure? (ii) What is the 
physical origin of the periodic periodic bright spots in the 
external ring? And, (iii) is the a coherent phenomenon, 
or a result due to unbalanced transportation? In Ref. d, 
a charge separated transportation mechanism was pro- 
posed. It gave a satisfactory explanation to the forma- 
tion of the exciton ring and the dark region between the 
internal and external rings. However, the remaining two 
questions are still not clarified satisfactorily. 

It is commonly believed that the periodic bright spots 
in the external ring are formed due to certain kind of in- 
stability. It generates the density modulation by a pos- 
itive feedback. Regarding the origin of the instability, 
Levite et al. considered that exciton states are highly 
degenerate and the instability comes from the stimulated 
scattering.^ When a local fluctuation occurs in the exci- 
ton density, it will lead to an increase in the stimulated 
electron-hole binding rate. The depletion of local carrier 
concentration then causes neighboring carriers to stream 
towards the point of fluctuation. On the other hand, 
Sugakov suggests that the instability is due to an attrac- 
tive interaction between the high-density excitonsi^ The 
creation of different structures, for example, the islands 
or the rings of the condensed phase, occurs due to the 
nonequilibrium state of the system connected with the 
flnite value of the exciton lifetime and the presence of 
pumping. Therefore, the appearance of the patterns is 
likely a consequence of self-organization processes in a 
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nonequilibrium systemiS 

To understand the subtle properties of indirect exciton, 
further experiments have been done and reported. For 
example, a novel method was proposed to demonstrate 
that cold exciton gases can be trapped by the laser in- 
duced trapping, which is similar to the trapping of BEC 
ultracold atoms.— More recently, an improved trapping 
technique was used to lower the effective temperature of 
indirect excitons.^^ With the even low temperature, four 
narrow PL lines have been observed, which corresponds 
to the emission of individual states of indirect excitons 
in a disorder potential. The homogeneous line broaden- 
ing increases with density and dominates the linewidth 
at high densities. 

In a recent paper (Ref. ITT| ). for the first time, Yang 
et al. measured the exciton PL energy along the cir- 
cumference of the ring. The most interesting result is 
that the average exciton energy depends on temperature 
nonmonotonically. With reducing temperatures, average 
exciton energy of the indirect exciton is lowered until the 
transition temperature, Ttr ~ 4 K, is met. Below Ttr^ 
the macroscopically ordered exciton state is formed but 
the average exciton energy in the ring actually increases 
with temperature being further lowered. In particular, 
the largest energy of single exciton is found to locate in 
the brightest regions. Due to these observations, Yang 
et al}^ argued that the interaction between exciton is 
repulsive. A numerical calculation seems to support this 
scenario. 

Another important experiment is the measurement of 
the coherence length. A Mach-Zehnder interferometer 
with spatial and spectral resolution was used to probe 
the spontaneous coherence in cold exciton gases, which 
are implemented experimentally in the ring of indirect 
exciton in CQWJ^ A strong enhancement of the exciton 
coherence length is observed at temperatures below a few 
Kelvin. The increase of the coherence length is correlated 
with the macroscopic spatial ordering of excitons. The 
coherence length at the lowest temperature corresponds 
to a very narrow spread of the exciton momentum distri- 
bution, much smaller than that of a classical exciton gas. 
It also shows that the apparent coherence length is well 
approximated by the quadratic sum of the actual exciton 
coherence length and the diffraction correction given by 
the conventional Abbe limit divided by ttj^ 

Whether the interaction between excitons in the 
macroscopically ordered state is attractive or repulsive 
and to what extent the interaction affects the forma- 
tion of the macroscopically ordered state remain an open 
question. The current paper attempts to clarify the 
above issues. In fact, experimental data had revealed sev- 
eral important clues. After analyzing the experimental 
facts, it is found that attractive interaction is dominant 
in indirect excitons and their nonequibirium distribution 
plays an important role in the macroscopically ordered 
state. Moreover, the uncertainty principle is shown to 
lead to a good account for the experimental phenomena 
and numerical simulations have confirmed that. 



This paper is organized as follows. In Sec.|TIl based on 
the uncertainty principle, we present a qualitative analy- 
sis of the temperature dependence of the average exciton 
energy in the external exciton ring. No interaction is 
considered in this section. In Sec. IIIIl the uncertainty 
principle is used again to explain the particle number 
density dependence of the PL FWHM (full width at half 
maximum) broadening and energy shift. The PL spectra, 
which are broadened first and then become sharper, pro- 
vide a strong evidence of the competition between a two- 
body attraction and a three-body repulsion. In Sec. IIVI 
we use a phenomenological nonlinear Schrodinger equa- 
tion, together with a temperature and energy dependent 
exciton distribution, to discuss the exciton energy spatial 
distribution in the external ring. The largest energy of 
single exciton is found to be in the brightest regions. In 
Sec. [Vj temperature dependence of the average exciton 
energy in the external ring is studied quantitatively using 
the approach proposed in Sec.|lVl The calculation shows 
that average exciton energy first decreases and then in- 
creases with increasing temperatures. In Sec. lVIl we give 
a demonstration on the exciton distribution in the laser 
induced trap case. Sec. I VIII is devoted to a brief sum- 
mary. 



II. TEMPERATURE DEPENDENCE OF 
ENERGY: QUALITATIVE ANALYSIS 

To understand the experimental phenomena, we first 
analyze the electron and hole creation, transportation, 
and the exciton formation qualitatively. As pointed out 
m Refs. dill, when electrons and holes are excited by 
lasers, they are hot electrons and hot holes initially. Since 
the drift speed of hot electrons is larger than that of 
hot holes (electron has a smaller effective mass), elec- 
trons and holes are indeed charge-separated at this stage 
(no true exciton is formed). Because hot electrons and 
holes have a small recombination rate and optical inac- 
tive, they can travel a very long distance from the laser 
spot. After a long-distance travel, hot electrons and holes 
collide with the lattices and are eventually cooled down. 
Due to the neutrality of the CQW, negative charges will 
slow down and accumulate far away from the laser spot. 
As a consequence, spatially separated electrons and holes 
form excitons and become optical active at the boundary 
of the opposite charges where they recombine and show 
a sharp luminescence ring. This kind of charge separated 
transportation mechanism has been used to interpret the 
exciton ring formation.i^ii^ 

It should be emphasized that, however, after the long- 
distance transportation and cooling, cooled electrons and 
holes will meet in the region of the external ring in the 
experiment of Butov et al.^ . However, in the experiment 
of Lai et al.^, they meet in the impurity potential well. It 
is believed that only cooled electrons and holes can form 
excitons. At this stage, charges are not separated, but 
coupled or bound together. The long lifetime of the ex- 
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citons means the lower electron-hole recombination rate. 
Therefore the elemental particles in the external ring'^ 
and in the impurity potential well^ are excitons. We thus 
argue that it is the interaction between excitons, rather 
than the unbalanced transportation, which leads to the 
nonhomogeneous density distribution associated with the 
complex PL patterns. 

When the experimental temperature is low (T < 1.4 
K) , average translational kinetic energy and in-plane mo- 
mentum of excitons are also low. In such case, excitons 
are optically active and the space-dependent PL inten- 
sity is proportional to the exciton number distribution. 
Since exciton number distribution is equal to the prod- 
uct of the probability distribution and the total excition 
number N, if N remains unchanged, PL intensity is di- 
rectly proportional to the probability distribution. In 
the following discussion, we simply take exciton density 
distribution as the PL distribution. Besides, since parti- 
cles can only move within the CQW, their motions are 
essentially two-dimensional (2D). 

The periodic array of beads in the external rings is a 
low-temperature phenomenon. This means that average 
exciton translational kinetic energy is low enough and 
its wave-particle duality is important. Therefore exciton 
energy is governed by the uncertainty principle. The un- 
certainty principle predicts that Ap h/Ar, where Ap 
is the momentum spread, Ar is spatial uncertainty, and 
h is the Planck's constant. Since exciton's momentum 
can be approximated by p ^ Ap, the corresponding en- 
ergy E ~ ?i^/[2m(Ar)^]. With these in mind, one is able 
to consider the temperature dependence of the exciton 
energy distribution, such as those reported in Ref. [TTI . 
When temperature is lower than Ttr — 4K, macroscopic 
ordered state forms. Exciton are confined into each bead 
of the external ring. Their spatial uncertainty Ar can be 
approximated by the diameter R of the bead (Ar ~ R). 
Experiments have observed that the lower the tempera- 
ture is, the higher contrast the pattern is and the smaller 
the R isr^ This implies that when temperature is lower, 
momentum uncertainty will be larger and hence the exci- 
ton energy will be higher. This gives a qualitative inter- 
pretation why the average exciton energy increases with 
decreasing temperatures. 

When temperature is increased to be higher than Ttr, 
macroscopically ordered state breaks down and the con- 
fined regions enlarge essentially. The uncertainty related 
energy will decrease in principle. However, when the 
temperature is increased further, the system undergoes a 
transition into the classical limit to which wave character 
of exciton is no longer important. In this case, the av- 
erage exciton translational kinetic energy is directly pro- 
portional to the temperature. Therefore, non-monotonic 
temperature dependence of the average exciton energy is 
closely related to the pattern formation, and has little 
to do with the interaction between excitons. In Fig. 3 
of Ref. [Til linear temperature dependence of the aver- 
age exciton energy is reported at T > Ttr- Their results 
confirmed that excitons behave classically in this limit 



to which E ~ ksT at 2D. Nevertheless, the slope of the 
linear temperature dependence was measured in Ref. [ll| 
to be about 1.4 x 10^'* eV k~^, slightly larger than the 
Bohzman constant, ks = 8.62 x 10"^ eV k'^. The dif- 
ference may be due to some effect of the interactions and 
extra degrees of freedom in addition to translational mo- 
tions. 



III. EXCITON PL SPECTRA: QUALITATIVE 
ANALYSIS 

For a long time, indirect exciton PL FWHM energy 
and the shift with increasing density in the presence of 
a random potential remain to be fully understood i^'^i^^ 
In a previous paper—, we have proposed a mechanism to 
study the macroscopically ordered exciton states. The 
most important point lies in the interaction which in- 
volves the competition between a two-body attraction 
and a three-body repulsion. Here we elaborate this point, 
along with the uncertainty principle, to illustrate the PL 
spectra. 

It is instructive to emphasize several important fea- 
tures in connection with our model. The idea first came 
from the density dependence of inhomogeneous exciton 
distribution. It is clear that the interaction between the 
indirect exciton is neither purely attractive, nor purely 
repulsive. If the interaction between the indirect exci- 
ton is purely repulsive, it will drive the exciton towards 
homogeneous distribution and the exciton cloud will ex- 
pand with the increase of the exciton number. On the 
contrary, if the interaction is purely attractive, the sys- 
tem is expected to collapse when the exciton density is 
greater than a critical value to which there is not enough 
kinetic energy to stabilize the exciton cloud. In addition, 
the case of a purely repulsive or a purely attractive inter- 
action is incapable to understand the experimental fact 
that the exciton cloud contracts first and expands later 
when the laser power is increased.^ 

The existence of the attractive interaction does not 
mean that the exciton state is unstable against the for- 
mation of metallic electron-hole droplet. The reason is 
that the repulsive interaction may dominate over the at- 
traction in that regime and keep the system stable. Many 
effects may contribute to the exciton interaction. The 
most important one is the exciton dipole-dipole interac- 
tion. Exciton behaves like a dipole, so a strong repul- 
sion will govern exciton interaction when two dipoles are 
aligned parallel. However, if the direction of two dipoles 
changes from aligned parallel to inclined, the attraction 
between the election of one exciton and the hole of the 
other exciton will dominate. This case can happen eas- 
ily when exciton density is low. In case of high densi- 
ties, due to the strong Coulomb interaction, the dipoles 
will tend to align parallel and consequently a repulsive 
interaction dominates. The other important interaction 
originates from the exchange effect. When two indirect 
exciton approach to each other, the exchange interaction 
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between two electrons, as well as the one between two 
holes become hxiportant. This may be another source of 
the attractive interaction between excitons. 

In fact, the complex exciton interaction which results 
from the competition of multi-effects, may be well de- 
scribed by the van der Waals form. It has pointed out 
that the effective interaction between exciton will be 
attractive when the separation between two exciton is 
about 3 to 6 exciton radii*^ In the current experiment, 
the exciton density is about 10^°/ cm^. For this den- 
sity, the average distance between the indirect exciton is 
about 100 nm. As the exciton Bohr radius qb is about 10 
~ 50 nm^, the average distance between excitons is about 
2 ^ 10 exciton radii. Thus it is reasonable to assume that 
the two-body interaction is in the attractive regime. In 
fact, the attractive interaction between the exciton has 
been considered as a possible candidate to describe the 
pattern formation observed by experimentSj^iiS 

In a previous paper—, taking into account the two- 
body attraction and the three-body repulsion interac- 
tions, we have proposed that in the dilute limit the be- 
havior of the exciton can be described by a nonlinear 
Schrodinger equationii 

^2 

- ^^'V^J + C^e- - 9in + 92n^)^j - E^^j, (1) 

where ipj and Ej are the j-th cigcnstate and eigenvalue, 
respectively. Vex is an external potential, and gi and 
g2 are (positive) coupling constants associated with two- 
body and three-body interactions, n — n{r) is the local 
density of exciton and at the mean field level, 

n{r) ^Y.'^,{E,M,{r)\^ (2) 
i=i 

where Af denotes the total number of bound states and r]j 
is the corresponding probability function associated with 
the energy level Ej . rjj satisfies the normalized condition 
= 1- Readers can refer to Ref. [13 for a 
detailed description of the model. 

In the present model, when macroscopically ordered 
state forms in the ring or in the impurity (disorder) po- 
tential well, exciton will distribute over the discrete en- 
ergy levels. In fact, four sharp peaks in PL spectra corre- 
sponding to the emissions of indirect exciton states have 
been observed recently in a so-called "elevated trap".^° 
We attribute the individual localized states as the con- 
finement due to both the elevated trap potential and the 
self-trapped potential (arising from the attractive inter- 
action). 

The present model can give a natural and self- 
consistent explanation of the particle density dependence 
of FWHM broadening and PL energy shift. For low den- 
sity exciton, the interaction between exciton is dominated 
by the attraction. With the increase of the particle den- 
sity (by increasing the laser power) , the low temperature 
exciton cloud will contract irrespective of localized or de- 
localized states. This means that the attractive interac- 
tion hampers the exciton motion. At the mean field level. 



this corresponds to the increase of the self-trapped po- 
tential. The larger the particle density is, the larger the 
attractive interaction is, and the exciton cloud will con- 
tract further. At low temperature, T < Ttr, uncertainty 
principle governs the exciton motion. Consequently, it 
results in higher mean exciton energy and larger energy 
dispersion. Higher mean exciton energy indicates blue- 
shift of PL peak and larger energy dispersion means that 
PL peak is more broadened (FWHM becomes large). 

When the particle density is further increased, three- 
body repulsion becomes more important and self-trapped 
interaction becomes weaker. The PL peaks will red shift 
and become sharp again. In fact, with the increase of 
the laser power in the low particle density regime, the 
phenomenon that PL spectra broaden first and then be- 
come sharper has been observed (see Fig. 3(c) of Ref. [Toi) . 
This, in addition to the support given by the experimen- 
tal data in Ref. 0, provides a further evidence that the 
interaction between exciton is likely the combination of 
a two-body attraction and a three-body repulsion. With 
increasing the laser power, blue shift of the PL peaks has 
been reported in Ref. [l^- However, the prediction of 
blue shift first and then red shift of the PL peak with the 
increase of the laser power has not been observed yet! 
Further experiments are in demand. 



IV. SPATIAL ENERGY DISTRIBUTION 

In order to study the spatial dependence of the en- 
ergy distribution in the external exciton ring, it is nec- 
essary to solve the nonlinear Schrodinger equation ([T]) 
numerically. In addition to the assumption of interac- 
tions, another important aspect lying in our model is 
the disequilibrium energy disequilibrium energy distribu- 
tions of indirect excitons. In fact, the indirect excitons in 
CQW involves both the (complex) energy relaxation and 
recombination processes. These have been studied by 
several experimental and theoretical groupsi^ It shows 
that for a relaxation process, when the exciton density is 
low {n <^ ttB is Bohr radius), the effects due to the 
exciton-exciton interaction and the exciton-carrier scat- 
tering can be neglected. In this case, the relaxation time 
is mainly determined by the scattering of excitons with 
acoustic phononsj^ In particular, at low bath tempera- 
tures {Tb < 1 K), this kind of relaxation rate decreases 
dramatically due to the so-called "phonon bottleneck" 
effects^ 

For the recombination process, on the other hand, be- 
cause exciton in the lowest self-trapped level are quantum 
degenerate, they are dominated by the stimulated scat- 
tering when the occupation number is more than a criti- 
cal value. Strong enhancement of the exciton scattering 
rate has been observed in the resonantly excited time- 
resolved PL experiments^. Therefore, even though the 
phonon scattering rate is still larger than the radiative 
recombination rate, thermal equilibrium of the system 
may not be reached. Essentially the distribution rij{Ej) 
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may deviate significantly from the usual Boltzmann or 
Bose like^^ 

A recent paper (Ref. [l3) reported that the observed 
narrow PL lines of indirect excitons in a disorder po- 
tential correspond to the emission of individual states. 
It has been shown that the intensity of PL spectra is 
roughly in proportion to the exciton number. When ex- 
citons distribute over some given energy levels, the PL 
spectra will exhibit strong peaks corresponding to these 
levels. Therefore the relative height of the PL peaks can 
be interpreted as the exciton distribution probability as- 
sociated with the energy level, rij{Ej). 

At the high temperature of 10 K, the PL peak inten- 
sity of the high excited state is higher than that of the 
ground state (see Fig. 2(c) of Ref. 10). This indicates 
that exciton have a larger probability distribution in the 
high excited state. With reducing the temperature to 
5.1 K, the PL peak intensity becomes lower for the high 
excited state while it becomes higher for the other three 
peaks associated with the two ground states. It seems 
that excitons at higher states will relax to lower states 
when the temperature is reduced. At the temperature 
lower than 5.1 K, the intensities are found to be almost 
identical for the four peaks. It indicates that excitons dis- 
tribute over the four discrete levels with almost the same 
probability. Furthermore, at even lower temperature of 
2.7 K, the low-energy peaks become very dim. This im- 
plies that phonon bottleneck effect is in effect for ground 
state and low excited states and in this regime, excitons 
in high excited states are having difhculty relaxing to the 
ground states. 

From the above analysis, one can generalize ri{Ej) — 
r]{Ej,T) to have a better description for the tempera- 
ture dependence of particle number distribution of bound 
states. As mentioned above, experiment seems suggest- 
ing that the weight of r]{Ej,T) will transfer from high- 
energy states to low-energy states when T is decreased. 
It also suggests that, for a given T, the weight of r]{Ej , T) 
is larger in higher-energy states than that in lower-energy 
states. With the above in mind, two other factors are also 
important in determining the actual form of rj{Ej,T). 
Firstly, exciton are Bose quasiparticle and at not too low 
temperatures, the distribution can be approximated by 
the (classical) Boltzman function, exp{—Ej/T), where T 
is considered as an effective temperature related to lattice 
temperature. The other important factor is the energy 
dependence of exciton luminance efficiency. Low-energy 
excitons is more optical active than that of high energy 
exciton. So the low energy excitons have a high lumi- 
nance efficiency which is contrary to the high-energy ex- 
citons. Thus low-energy state will have relatively smaller 
exciton distribution taking into account the luminance 
efficiency. One can use exp(i?j/i?o) to describe qualita- 
tively this effect, where Eq is an effective energy scale 
related to the exciton life time. As a consequence, tem- 
perature and energy dependent exciton distribution is an 
outcome of the competition between the above two ef- 
fects. We then take the following phenomenological form 
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FIG. 1: Energy dependent exciton distribution. Solid line is 
obtained from the phenomenological formula ((3]) with chemi- 
cal potential fi = 1.56 eV and a = 0.025 eV~^. Circle line cor- 
responds to the exciton PL spectra (in arbitrary unit) taken 
from Fig. 2(a) of ^ld\. For comparison, a constant is multi- 
plied to the solid line. 

for the weight 

viEj,T) = Ce^p[-iE,-^l)/T]el,p[iE,-^i)/Eo] 
= Cexp[a{E,~fi)], (3) 

where C is the normalization factor, ^ is the chemical 
potential, and a = I/Eq — 1/T. Typically Eq < T and 
hence a > 0. When temperature and laser power re- 
main unchanged, fi and a will remain constantly. While 
increasing the temperature, a will increase correspond- 
ingly. Thus a can be considered as an "effective temper- 
ature" . 

Equation ([3]) basically gives an appropriate temper- 
ature and energy dependent probability density distri- 
bution for indirect exciton. In Fig. [l] we compare the 
theoretical curve [based on Eq. (15] )] with the experimen- 
tal PL spectra reported in Ref. [lO. The experimental 
data (circles) in Fig. [T] are in arbitrary unit. For com- 
parison, we have multiplied a constant to the probability 
density. It should be mentioned that we have tried var- 
ious distribution functions and found that Eg. ([3]) gives 
the best fitting. In our earlier study (Ref. 17), we have 
assumed that exciton are distributed at different energy 
levels with almost the same probability. Strictly speak- 
ing, it corresponds only to the low temperature limit in 
Eq. @. 

Numerically it is convenient to do the following scal- 
ing: 2Pj{r)/VN iPj{r), Ngi gi, and N^g2 ^ 52, 
such that Eq. ([T]) remains the same form. In this case, 
n(r) becomes the probability density which satisfies the 
normalization condition Jgn{r)dS = 1. After further 
rescaling ^j{r)aph — >■ i/'j(r) and r/(TpL r, Eq. ^ is 
reduced to 

- ^^^V'j + (vex - ain + a2n^)tpj = e-jtpj, (4) 

where v^x = V^x/e, ai = 51/ (c^pl<^)' = 52/ {<^pl^), 
and Ej = Ej/e. Here e = h^/m*apj^ with trpL being the 



root-mean-square radius of the exciton cloud observed by 
photoluminescence. With the above scaling, it is found 
that a\/a2 = 91/92^ — 5i'ti*(Tpl/52^^, which is a con- 
stant for an explicit sample. 

Taking into account the experimental facts, '^'-1 two im- 
portant points should be clarified, (i) The exciton pat- 
terns are fully determined by its self-trapped interactions. 
External potential V^x is not the main cause for complex 
exciton patterns. Thus we set v^x = for simplicity, 
(ii) When an electron and a hole form an exciton at low 
temperature, it is believed that their kinetic energy is 
low and all exciton are and all exciton are self-trapped. 
Particles with energy greater potential energy can not be 
bounded in the self-trapped well and should be ruled out 
in the calculations. 

The mean kinetic energy of excitons can be given by 



Ek= j J dxdy Ek{x,y), 



(5) 



where Ek{x,y) is the mean local (space-dependent) ki- 
netic energy density. Considering the probability func- 
tion of each level, rjj = ri{Ej,T), Ek can also be given 
by 



Ek 



M 

^rijEkj, 



(6) 



where Ekj is the kinetic energy associated with level j. 
In fact, 



Ek 



kj 



dxdy ilj*{x,y)V'^ipj{x,y) 



dxdy 



d'<Pjix,y) 



dx 



diljj{x,y) 



dy 



,(7) 



with ipjlxjy) the wave function of level j. In obtain- 
ing the 2nd line of Eq. (O , an integration by parts and a 
boundary condition that wave function vanishes at the in- 
finity (bound states) are applied. By comparing Eqs. ([5])- 
([7|), the mean spatial energy density is obtained to be 



A/" 



Ek{x,y) = ^?yj 



d'ipj{x,y) 



dx 



dtpj{x,y) 



dy 



(8) 



Besides, one can also obtain the angle dependent kinetic 
energy density via 



Ek [9) 



dxdy Ek{x, y) 



(9) 



{x,y)<£{eM+Ae} 



with tan 9 = y/x. In Fig. [D^a), we show both the par- 
ticle density distribution, n(r), in the exciton ring and 
the angle-dependent kinetic energy distribution, Ek (9). 
Obviously the maxima of Ek (9) are located at the center 
of the beads. Fig. [2Kb) shows the spatial dependence of 
the kinetic energy distribution, Ek{x,y). One interest- 
ing feature is that Ek{x,y) is zero at the center of the 
beads. This can be elucidated by future experiments of 
finer resolution. 




FIG. 2: (a) Formation of the macroscopically ordered state of 
the exciton ring at very low temperatures. The corresponding 
angle-dependent energy distribution Ek{9) (solid line) is also 
shown, (b) Spatial distribution of the kinetic energy, Ek{x, y). 
Parameters used are ai — 250, 02 — O.OOlaf, effective tem- 
perature a = 0.002, and chemical potential fj, = 0. 



V. TEMPERATURE DEPENDENCE OF 
ENERGY: QUANTITATIVE RESULTS 

In section |IT1 with the uncertainty principle, a qualita- 
tive analysis was made on the temperature dependence of 
exciton average energy in the external ring. Here we solve 
the nonlinear Schrodinger equation ([l} , together with the 
disequilibrium distribution rij{Ej,T), to study the tem- 
perature dependence quantitatively. As mentioned be- 
fore, a can be treated as an effective temperature in the 
present case. Temperature (a;)-dependent exciton ring 
patterns are shown in Fig.[3Ka)-(c), while Fig.[2Kd) shows 
the nonmonotonic temperature dependence of the total 
energy, Ek ^ ! Ek{9)d9. 

It should be noted that our theory applies only to low 
temperatures when excitons are in or near the condensed 
state. When temperature is above the critical point, such 
as Ttr — 4 K presented in Ref. [13, excitons will trans- 
fer from condensed to non-condensed phase. Above Ttr 
in the non-condensed phase, the exciton liquid may be 
"boiling" and its spatial distribution could change with 
the time. At such high temepratures, our theory will 
break down to which numerical calculation turns no con- 
vergent solutions. In fact, the exciton pattern shown in 
Fig. [SKc) is just a snapshot of time-dependent exciton 
number distribution. When a > 0.08 in our numerical 
simulation, the exciton behavior can be well described 
by the classical statistical physics and a linear formula 
£ = ksT is plotted in Fig.[3jd) for illustration purpose. 

A brief summary is in order here. Initially, if exciton 
are uniformly distributed on the external ring, attractive 
interaction will drive the exciton to approach each oth- 
ers. When the local density reaches a critical value, the 
kinetic energies of the exciton will drive the high-density 
excitons to diffuse. At the same time, the repulsive in- 




FIG. 3: Formation of the macroscopically ordered state of the 
exciton ring at different effective temperatures: (a) a = 0.002, 
(b) a = 0.004, and (c) a = 0.008. (d) Temperature (a) 
dependence of the total energy Ek- All other parameters are 
taken to be the same as those in Fig. [S] 

teraction, which increases with the density, will hinder 
further increase of the exciton density. As a result, an 
array of clusters on the external ring forms. The size of 
these clusters is thus determined by three factors, i.e., the 
attractive interaction, the three-body repulsive interac- 
tions, and the kinetic energies of the excitons associated 
with the temperatures. This is the basic idea behind our 
model on the macroscopically ordered states of excitons. 

Finally some remarks on the temperature effect are 
given. When the bath temperature is low, cooled ex- 
citon have relatively low momenta and the self-trapped 
interaction is able to confine most of the excitons. How- 
ever, in the low momentum case, cooling efficiency is low 
while luminous efficiency is high, excitons can not reach 
the thermal equilibrium state. At the meantime, due to 
the competition between the self-trapped and the kinetic 
energies, complex exciton patterns occur (as discussed 
above). With increasing the temperature, the exciton 
can not be fully cooled and correspondingly self-trapped 
interaction confines only part of the excitons. The at- 
tractive interaction also cannot compensate the exciton 
kinetic energy and excitons will distribute homogenously 
in a 2D plane. In this case, the pattern is washed out. 
If the temperature is higher than the indirect exciton 
binding energy ^ 3.5 meV^^i^^ most of excitons become 
ionized and are in a plasma state. No pattern can be 
observed in this case. 



VI. NUMBER DISTRIBUTION WITH LASER 
INDUCED TRAP 

Similar to the studies of ultracold alkali atoms and 
molecules, the laser-induced trapping was proposed 
and demonstrated for a highly degenerate Bose gas of 
excitonJ?. An important advantage of laser trapping is 
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(a) _ (b) m 




FIG. 4: The exciton distribution in laser-induced trap for 
(a) low (ai = 15), (b) intermediate (ai — 23), and (c) high 
(fli = 35) densities. The other parameters are taken to be 
same as those in Fig. m 

that it is possible to control the trap in situ by varying 
the laser intensity in space and time. A commonly seen 
phenomenon for highly degenerate excitons is its annual 
particle distribution in the laser-induced trap. We argue 
that this kind of distribution is similar to that found in 
an impurity potential reported in Ref. 0. In our pre- 
vious work, a detailed discussions was already given to 
this kind of trapped exciton distribution?!? Since laser 
induced trapping can confine higher density exciton, it is 
hoped that the experiment can also be done at very high 
densities to explore the many-body physics of excitons. 

Using the temperature and energy dependent distri- 
bution given in Eq. ([3]) together with a two-dimensional 
pseudopotential 

/ \ / -5 for r < o-pL 
^-W = \ otherwise, ^^O) 

we solve the nonlinear Schrodinger equation ([!]) for exci- 
ton number distribution ranging from low to high densi- 
ties. The results are presented in Fig. SI As mentioned 
before, (Tpl is the root-mean-square radius of the exciton 
PL pattern. To compare to experiments, we take ai = 
15, 23, and 35 with 02 = 0.004af accordingly. As shown 
in Fig. m the excitons distribute annular ly at low densi- 
ties and at high densities, some fine structures develop at 
the center of the ring, which is consistent with the experi- 
ments. Readers can refer to Ref. [13 for the reason behind 
the annular distribution and some related physics. 

VII. SUMMARY 

In summary, the uncertainty principle is used to ana- 
lyze the spatial and temperature dependence of the av- 
erage exciton energy distribution in a macroscopically 
ordered state of the exciton. the competition between 
a two-body attraction, a three-body repulsion, and the 
kinetic energy plays a crucial in determining the behav- 
iors of exciton distributions. Numerical simulation of the 
corresponding nonlinear Schrodinger equation seems to 
confirm the analysis. Nevertheless, the reason of form- 
ing macroscopically ordered exciton states may be more 
complex than what is expected. Full understanding of 
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the exciton interactions requires a reliable many-body 
calculation beyond the mean-field approximation, which 
is obviously not feasible at the moment. In order to re- 
alize the exciton BEC phenomenon, it is in demand to 
obtain the exciton with even lower combination rate and 
shorter relaxation time. 
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